Magnetic resonance imaging apparatus and image processing method

ABSTRACT

A magnetic resonance imaging apparatus having control means for fetching a plurality of image data and obtaining an image by a calculation. The control means performs an unwrap processing to a plurality of regions for removing a main value of a phase rotation amount so as to obtain a distribution diagram of a phase rotation amount caused by magnetostatic unevenness, obtains a plurality of evaluation values based on the phase value between a plurality of regions, selects a desired evaluation values based on the phase value from the plurality of evaluation values, and adjusts the phase values between the regions in accordance with the desired evaluation value.

TECHNICAL FIELD

[0001] The present invention relates to a magnetic resonance imaging apparatus (hereinafter referred to as MRI apparatus); and, more particularly, the invention relates to a technique for accurately obtaining a water/fat separated image of each of two or more island-like regions included in one image.

BACKGROUND OF THE INVENTION

[0002] When a magnetic resonance imaging apparatus is used to obtain an image, images showing the contrast between various tissues can be obtained by changing parameters of the time of echo (TE), the time of repetition (TR), and performing image calculations. The methods of obtaining an image having the contrast include the Dixon method. The Dixon method is a method of obtaining a water/fat separated image by utilizing the phase change of water MR signals and fat MR signals in one region.

[0003] Recently, in clinical settings, opportunities of imaging both right and left crural regions at the same time, or of imaging a transverse cross section of brain-basilar region of head, or of imaging a coronal cross-section of abdomen are increasing in diagnosing from the water/fat separated image. A conventional technique of accurately obtaining a water/fat separated image when the imaging region is separated into a plurality of island-like portions since an object to be examined includes a plurality of spatially-separated portions, or an object includes portions without signals, such as air is disclosed in “A New Two-Dimensional Phase Unwrapping Algorithm for MRI Images”; M. Hedley and D. Rosenfeld; Magnetic Resonance in medicine, Vol.24, 177-181 (1992). According to this technique, when there are two island-like regions, only the phases of the points of the first region and the second region closest to each other are simply compared. However, since, with the above technique, it is difficult to obtain a desired image of a portion where there is a local rapid change of phases, or where the island-like regions are greatly separated, the water image and the fat image cannot be separated, and portions where the water signals are extracted and portions where the fat signals are extracted are mixed up in one image.

SUMMARY OF THE INVENTION

[0004] The object of the invention is to provide a magnetic resonance imaging apparatus that can suppress the generation of artifacts between imaging regions even when one image includes two or more regions that are separated from each other. According to this invention, a magnetic resonance imaging apparatus is disclosed, which includes control means for acquiring data of a plurality of images having various echo times, and obtaining an image by calculation, wherein the control means adjusts phase values between the respective portions using an evaluation function based on the phase values of the respective portions.

BRIEF DESCRIPTION OF THE DRAWINGS

[0005]FIG. 1 is a block diagram showing the structure of an MRI apparatus;

[0006]FIG. 2 is a diagram illustrating the conventional static magnetic field inhomogeneous distribution within an FOV;

[0007]FIG. 3 is a flow chart illustrating the cross-region unwrapping algorithm according to the invention;

[0008]FIG. 4 is a diagram showing an image of an object to be imaged, which includes a plurality of portions that are spatially separated;

[0009]FIG. 5 is a graph showing the phase profile of an image of an object to be examined, which includes a plurality of portions that are spatially separated;

[0010]FIG. 6 is a graph illustrating a static magnetic field inhomogeneous distribution within the FOV of an object to be examined, which includes a plurality of portions that are spatially separated; and

[0011]FIG. 7 is a graph comparing evaluation values according to the invention.

[0012] The numerical reference 101 represents an object to be examined; the numerical reference 102 represents a static magnetic field magnet; the numerical reference 103 represents a gradient magnetic field coil; the numerical reference 104 represents an RF coil; the numerical reference 105 represents an RF probe; the numerical reference 106 represents a signal detection unit; the numerical reference 107 represents a signal processing unit; the numerical reference 108 represents a display unit; the numerical reference 109 represents a gradient magnet field power source; the numerical reference 110 represents an RF transmission unit; the numerical reference 111 represents a control unit; and the numerical reference 112 represents a bed.

BEST MODE FOR CARRYING OUT THE INVENTION

[0013] The structure of an MRI apparatus will be described with reference to FIG. 1. The apparatus includes a magnet 102 for generating a static magnetic field around an object 102 to be examined; a gradient magnetic field coil 103 for generating gradient magnetic fields in the space; an RF coil 104 for generating an RF magnetic field in the space; and an RF probe 105 for detecting MR signals generated from the object 101. The gradient magnetic field coil 103 consists of three gradient magnetic field coils, respectively in the X, Y, and Z directions, and each of them generates a gradient magnetic field in response to signals from a gradient magnetic field power source 109. The RF coil 104 generates an RF magnetic field in response to signals from an RF transmission unit 110. Signals from the RF probe 105 are detected by a signal detection unit 106, subject to signal processing in a signal processing unit 107, and are converted into image signals by calculation. The image is displayed on a display unit 108. The gradient magnetic field power source 109, the RF transmission unit 110, and the signal detection unit 106 are controlled by a control unit 111. The time chart of the control is generally referred to as a pulse sequence. A bed 112 is designed for the object to be laid on.

[0014]FIG. 2(b) is a static magnetic field inhomogeneity map in imaging an object 208, and FIG. 2(a) is a graph where phase values along the line A-A′ are plotted. In FIG. 2(a), the horizontal axis represents the position, and the vertical axis represents a rotation amount of the phases caused by the inhomogeneity of the static magnetic field. The graph shows a distribution 201 of the phase rotation amount caused by the static magnetic field inhomogeneity inside the object. Since the signal component outside the object 208 is substantially 0, and the phase value becomes measurement noise (random phase), a threshold absolute value of the signal is provided. The phase values of a portion where the signal value is smaller than a constant value are not used.

[0015] In FIG. 2(a), phase values 202 that are +π or more are wrapped to have a value 203. Values 204 and 205 that are −π or less are also wrapped to have values 206 and 207 between −π and +π. Thus, the phase value discontinuously jumps from −π to +π (the broken lines in the figure) in the wrapped area. However, the actual magnetic field does not jump discontinuously in such a manner. Consequently, the unwrapping processing mentioned in U.S. Pat. No. 5,701,074 is performed so as to obtain a smooth static magnetic field inhomogeneity map, eliminating the discontinuous jump.

[0016] Two separated objects to be unwrapped are shown in FIG. 6. FIG. 6(b) shows a static inhomogeneous magnetic field map, where the region having the larger area is a first region 610 and the region haivng the smaller area is a second region 611. FIG. 6(a) is a diagram where the phase values along the line A-A′ in FIG. 6(b) are plotted. In FIG. 6(a), the horizontal axis represents the position, and the vertical axis represents an amount of phase rotation caused by the static magnetic field inhomogeneity. Between the first region 610 and the second region 611 exists an air gap 613.

[0017] The cross-region unwrapping algorithm, shown in FIG. 3, is applied to such a processing object. Fitting is performed for calculating an approximate function of the phase values of the first region, the phase values of the second region is calculated, and then the unwrapping is performed.

[0018] The cross-region unwrapping algorithm, shown in FIG. 3, is also applied. First, the first region 610 is subject to the cross-region unwrapping (310). Next, two-dimensional fitting is performed only on the phase values of the first region 610 after the unwrapping (301). According to this embodiment, the fitting degree is four, for example, and an evaluation value f0 is calculated from an evaluation function based on the fitting result in the second region 611 (303). The evaluation function is a square-sum (SD: Standard Deviation) of a difference between the above fitting result and the phase values of the second region before the fitting, for example.

[0019] Then, the cross-region unwrapping is performed on the second region 611 (304). Next, the fitting is performed on the phase values of the first region 610 after the unwrapping (301) and the phase values of the second region after the unwrapping (304). Then, an evaluation value f1 is calculated from the evaluation function in the same manner (306). In the step 306, the evaluation value is calculated using the square-sum of the difference between the phase values before and after the fitting, as in the step 303.

[0020] Then, the phase values of the second region after the cross-region unwrapping is given an offset of +2π (307), the fitting is performed in the same manner (308), and thus, the evaluation value f2 is calculated (309). Further, the phase values of the second region after the cross-region unwrapping is given an offset of −2π (310), the fitting is performed (311), and an evaluation value f3 is calculated in this state (312). Here, the thus-calculated evaluation values f0, f1, f2, and f3 are compared, and it is judged which value is most suitable (313).

[0021] In another method, the unwrapping is not performed firstly on the first region. The fitting is performed on the phase values of the first region and an evaluation function of the second region is calculated. Then, the unwrapping is performed on the phase values of the first region and the second region, and the evaluation functions of the respective regions are compared. Then, it is judged which value is most suitable. The manner of this comparison is the same as that in FIG. 3.

[0022] In the judgment, the smallest square-sum is selected, for example. This process will be described with reference to FIG. 4 and FIG. 5. FIG. 4 shows an image of an object to be examined having two portions 401 and 402 that are spatially separated, and FIG. 5 shows a profile of phases of this image along the line B-B′. Areas 403 and 404 shown in FIG. 4 correspond with areas 501 and 502 shown in FIG. 5. A curve 503 represents the result of the phase unwrapping in the area 501. The curve 504 represents the result of the cross-region unwrapping in the area 502 using the above-mentioned cross-region unwrapping processing. The curve 504 has an offset of 2nπ, which means that the actual phase is represented by the curve 505. A curve 506 and a curve 507 represent the result of the fitting using the curve 503 and the curve 504, and using the curve 503 and the curve 505, respectively.

[0023] Here, when the curve 504 and the curve 506, and the curve 505 and the curve 507 are respectively compared, the trajectory of the curve 505 substantially coincides with that of the curve 507. On the other hand, the trajectory of the curve 504 is greatly different in values from that of the curve 506 in the area 502. Accordingly, it is judged that the state, where the square-sum of the difference between the phases after the cross-region unwrapping and the fitting result is minimum, is most suitable.

[0024] As a result of judgment in the step 313, if it is judged that the evaluation value f0 is most suitable; that is, the phase values before the cross-region unwrapping in the second region are most suitable, the values are changed back to the values before cross-region unwrapping in the second region (314). If it is judged that the evaluation value f1 is most suitable; that is, the phase values after cross-region unwrapping in the second region are most suitable, the value after cross-region unwrapping in the second region is employed (315). If it is judged that the evaluation value f2 is optimum; that is, the phase values that are givin an offset of +2π after the unwrapping in the second region are most suitable, the evaluation value fn is calculated while changing the offset value to +4π, +6π, . . . , so as to find the optimum offset value, and the thus-calculated offset value is employed (316). Similarly, if the evaluation value f3 is most suitable; that is, the phase values that are givin an offset of −2π are most suitable, the values of the evaluation function are calculated while changing the offset value to −4π, −6π, . . . , so as to find the optimum offset value, and the thus-calculated offset value is employed (316).

[0025] The process of the above-described judgment is schematically shown in FIG. 7, wherein the graph (a) shows the result of the steps 301 and 302. When the difference between a heavy line (actual phase value) and a thin line (fitting result) in the zone 701 is large, the evaluation value f0 becomes large. The graph (b) shows the result of the steps 304 and 305, wherein the values of the heavy line are very close to those of the thin line 704, and the evaluation value f1 becomes small. The graphs (c) and (d) show the result of the steps 307 and 308, and the result of the steps 310 and 311, respectively. The difference between the heavy line and the thin line 706 in the zone 705 and that between the heavy line and the thin line 708 in the zone 707 are large, and the evaluation values f2 and f3 also become large. In the example shown in FIG. 7, the evaluation value f1 is minimum, and the step 315 is found most suitable as a result of the comparison and judgment in the step 313.

[0026] When the above-described processing was applied to simultaneous imaging of left and right crural regions (transverse cross section or coronal cross section), multi-slice imaging of an abdominal transverse cross section or an abdominal coronal cross section, imaging with a large FOV, imaging of a head tranverse cross section, imaging of ankles, and the like, and three-hundred or more images were taken, the cross-region unwrapping was normally performed on all images and water/fat separated images without artifacts were obtained. Further, this method can be applied to both the three-point Dixon method of the gradient echo type and the three-point Dixon method of the spin echo type.

[0027] Incidentally, another evaluation function and judging formula may include the following examples:

[0028] (a) It is judged whether, when the fitting is performed on the phase rotation amount distribution map based on the static magnetic field inhomogeneity using a polynomial equation, all coefficients of the respective highest degree terms are the predetermined threshold value or less or not. If they are the threshold value or less, they are regarded to be right.

[0029] (b) In the above-described fitting, the square-sum of the coefficients of the highest-degree terms is calculated, and the minimum value among them is employed.

[0030] Here, (a) is based on the reasoning that a static magnetic field generated by a magnet gently changes; a phase map thereof does not have local change; and a coefficient of a high-degree term does not have a large value in the fitting result. In the actual processing, quadric fitting is performed, wherein the coefficients of quadric terms (x4, x3y, x2y2,xy3, and y4) are regarded to be right if they are a certain value or less. As a result of experiments, the processing can be stabilized, to some extent, by using this evaluation function. Incidentally, the optimum threshold value changes depending on the FOV and the portion imaged.

[0031] (b) is an improvement of (a), wherein a threshold is not provided in advance, but values of the evaluation function in the respective states are compared so as to judge which state is most suitable. With this method, a stablar result than that brought by (a) can be obtained. In this case, for the same reason as (a), the state where the value of (b) is minimum was regarded to be most suitable. In the same manner as (a), the processing result depended on imaging conditions.

[0032] In the above-described embodiments, among the island-like regions, the designation of the first region and the second region may be determined by a method other than comparison of their size. For example, the island-like region which is nearer the center of the image may be made the first region, and the other island-like region in the peripheral part may be made the second region. That is, it is important for performing the unwrapping rapidly and stably to start the unwrapping from the region which has a larger area and is nearer the center of the image in the processing order.

[0033] To calculate the evaluation function, the square-sum of the difference between the phases before and after the fitting in the second region is calculated. However, it also may be preferable to calculate the square-sum of the difference between the phases before and after the fitting in both the first and the second regions. In this case, not only the local fitting condition but also the fitting effect over the whole object is evaluated. However, since the area of the evaluated portion increases and the calculation of the evaluation value takes more time, more rapid processing is necessary.

[0034] Further, according to the above-described embodiment, the number of island-like regions to be unwrapped is two; however, the embodiment also can be applied when three or more island-like regions exist. In this case, a stable result can be obtained in a short time by employing the method where, after the cross-region unwrapping between the first region and the second region is completed, the cross-region unwrapping is performed on a third region. The evaluation value for the third region is calculated by the cross-region unwrapping between the first region and the second region, that is, from the evaluation function based on the phase values that are adjusted between the first region and the second region. Further, when a fourth region exists, the evaluation value is calculated from the evaluation function based on the phase values adjusted among the first to the third regions. Further, it is also preferable that, in the same way as evaluating the fitting effect of the combined first and the second regions, a square-sum of the difference of the phases before and after the fitting in the combined first and the ath (a is an integer, which is three or more) regions, so as to evaluate the fitting effect in the whole body of the object.

[0035] It is necessary to fully consider how large the island-like region to be unwrapped can be, for stably operating the processing. This can be easily judged by considering how advantageous it is to include a small island-like region, that is located at an edge of the image, in the calculation of the static magnetic field map. Typically, it is desirable that, in an image of 256×256, an island-like region including 100 pixels or less is disregarded and not used for determining the static magnetic field map. The reason is that a very small region might be a bloodstream artifact or a body-movement artifact, the phase values of which are unlikely to reflect the static magnetic field distribution. Even if an evaluation value (square-sum) is calculated for a small island-like region, there is no guarantee that the calculated value is accurate. Further, since the effect brought by such a small region on the entire fitting function is slight, the small region does not have much effect on the result of the static magnetic field distribution. Further, when the number of the island-like regions increases, it is necessary to consider extension of the calculation time. Further, the location is a criterion because the clinical importance of objects near the edge of the image is less. It is also effective to limit the number of island-like regions to be unwrapped. For example, only five island-like regions in descending order of the number of pixels are chosen for calculation. Such limitation has the advantage of limiting the calculation time within a certain range.

[0036] Further, according to the above embodiment, the quadric fitting function is used; however, the degree of the function may be varied. In this case, the following needs to be considered in the determination of the degree: when the degree is small, the static magnetic field inhomogeneity map of the acquired fitting result forms a smooth line, and the processing result is stabilized. However, there is a disadvantage that, with a low degree, the function cannot follow the local fluctuation of the magnetic field. For example, it is well known that in an MRI apparatus employing an open magnet, the magnetic field uniformity is generally deteriorated dramatically at the position farther from the magnetic field center. Therefore, when the size of the field of view (FOV) is increased and imaging is performed at a position other than the magnetic field center, the fluctuation of the static magnetic field cannot be accurately calculated. As a result, water and fat are mixed in a water/fat separated image whose static magnetic field is corrected, whereby the accuracy of clinical diagnosis is dramatically deteriorated. Further, it is known that, in a region where air space is complicatedly confusing, such as a nasal cavity, or a paranasal cavity, the static magnetic field locally fluctuates. In such a region also, the fitting cannot be accurately performed with a small degree. On the other hand, although those problems are solved by increasing the fitting degree, the fitting result is diverged around the object and become unstable. Further, the calculating time of the fitting is likely to be extended. In light of this experience, the inventors obtained a good result using a quadrate fitting function. Concretely, in the open MRI apparatus, the processing can be stably performed until when the FOV is 320 mm.

[0037] The reason why island-like regions to be unwrapped appear is, as described above, that pixels having a small absolute value are regarded as noise and are excluded from the object calculated. Although such handling is is regarded as quite natural for people in the same business, it is necessary to fully consider where to set the threshold value for stably performing processing. Typically, 5% to 40% of the maximum value of the signal in an image is desirable. If the threshold value is too large, the area of each island-like region becomes small, and the accuracy of the static magnetic field inhomogeneity map to be calculated is deteriorated since the map is created with little data. Further, the number of the island-like regions is likely to increase. If the threshold is small, the number of regions to be separated is reduced; however, since noise components are mixed in, the accuracy of the static magnetic field inhomogeneity map is deteriorated. According to the investigation by the inventors, it is preferable to set the threshold value in the range from 10% to 30% for stably performing processing.

[0038] The number of island-like regions also depends on the characteristics of the echoes to be measured. When a phase map is created by calculation using in-phase echoes of water and fat, the number of the island-like regions is reduced. However, when out-phase echoes of water and fat is used in the calculation for creating the phase map, the number of the island-like regions increases. The reason is that, in pixels where water and fat mix in an out-phase image, the respective signal components cancel each other, and thus the signal value becomes small. Consequently, to stably create the static magnetic field inhomogeneity map, it is most preferable to use both in-phase echoes and the phase unwrapping between island-like regions according to the invention. The stability and accuracy of the processing also can be greatly improved in comparison with the conventional processing, even when the invention is used with out-phase echoes, as described above. Further, the invention also may be applied when a static magnetic field inhomogeneity map is calculated from the phase difference between in-phase echoes and out-phase echoes. As a result of actually applying the invention to such data, the unwrapping processing on an image of brain-basilar region, an abdominal transverse cross section, and an abdominal coronal cross section was dramatically improved.

[0039] It is also important in which order the island-like regions are unwrapped. Experimentally, it is preferable to start the unwrapping from a region near the center of the image. The reason thereof is that, since the image center often coincides with the magnetic field center, the static magnetic field uniformity is high and the number of unwrapping in air gap is small, so that the obtained result is likely to be accurate. Otherwise, the unwrapping is preferably performed from a larger island-like region. The reason is that, since the accuracy of the evaluation value is high in a large island-like region, as described above, judgment is performed accurately, and the obtained result is likely to be accurate. Further, the unwrapping is preferably performed from an island-like region that has a smaller air gap with a region that has already been unwrapped, since the number of the unwrapping operations in the air gap is correspondingly small, and the obtained result is likely to be accurate. That is, the unwrapping is applied, starting with the region where a more accurate result is likely to be obtained. Since all of those processing orders can be described with simple equations, the automation of the processing is not interrupted.

[0040] According to the above embodiments, the fitting is used at a preliminary stage of calculating the evaluation value. Although a method of comparing the fitting result with the original data is one of the appropriate methods, the invention is not limited thereto. That is, it is also preferable that, when the unwrapping is performed between island-like regions, a plurality of unwrapping conditions are hypothesized, judgment is made on the respective results based on these conditions using an evaluation function, and an evaluation value that is likely to be most suitable, that is, the most suitable unwrapping condition, is employed.

[0041] Since the present invention is constructed as described above, even when a region to be unwrapped is separated into island-like regions, the cross-region unwrapping is properly performed. Therefore, since the phase rotation amount is accurately calculated, artifacts can be suppressed so that a region where a water signal is extracted and a region where a fat signal is extracted do not mix in one image. 

What is claimed is:
 1. A magnetic resonance imaging apparatus comprising control means for acquiring data of a plurality of images having various echo times and obtaining an image by calculation, said control means adjusting phase values between a plurality of regions by using an evaluation function based on phase values in the regions.
 2. A magnetic resonance imaging apparatus according to claim 1, wherein said control means executes unwrapping processing for removing the principal value of a phase rotation amount on a plurality of regions and adjusts phase values between the regions.
 3. A magnetic resonance imaging apparatus according to claim 2, wherein said control means calculates a plurality of evaluation values based on phase values between the plurality of regions, selects a desired evaluation value from among the plurality of evaluation values, and adjusts phase values between the respective regions based on the desired evaluation value.
 4. A magnetic resonance imaging apparatus according to claim 3, wherein said control means calculates the evaluation values based on phase values that have been subject to fitting after the unwrapping a distribution map of the phase rotation amount and also on phase values before the fitting.
 5. A magnetic resonance imaging apparatus according to claim 3 or 4, wherein said control means selects a minimum evaluation value from among the plurality of evaluation values.
 6. A magnetic resonance imaging apparatus comprising: static magnetic field generating means for generating a static magnetic field in space where an object to be examined is placed; RF pulse applying means for repeatedly applying an RF pulse that generates nuclear magnetic resonance phenomenon to the object; gradient magnetic field generating means for applying gradient magnetic fields in the slicing direction, the phase-encoding direction, and the readout direction, respectively; receiving means for detecting a plurality of echo signals from the object; and control means for controlling each means so as to generate a plurality of echo signals having various echo times within one repetition time of the RF pulse, apply a readout gradient magnetic field in parallel with the generation of each echo signal, and vary the phase encoding at each repetition, wherein said control means executes first unwrapping processing on the object having separated regions for adjusting phase values in each region, and executes second unwrapping processing for adjusting phase values between each region, in which a plurality of patterns of phase value adjustment are created, the desirable pattern is selected from among the plurality of patterns, and the phase values are adjusted by the selected pattern.
 7. An image processing method for image-processing a plurality of images having various echo times so as to obtain an image where water and fat are separated, comprising steps of: (a) performing the unwrapping processing on an image including a plurality of separated regions for calculating a distribution map of a phase rotation amount caused by the static magnetic field inhomogeneity in each region; and (b) calculating a plurality of parameters for adjusting phase values between each region, selecting a desired parameter from among the plurality of parameters, and performing the unwrapping processing between each region based on the selected parameter.
 8. An image processing method for image-processing a plurality of images having various echo times so as to obtain an image where water and fat are separated, comprising steps of; (a) performing the unwrapping processing for calculating a distribution map of a phase rotation amount caused by the static magnetic field inhomogeneity in a region to be standard on an image including a plurality of separated regions; (b) performing fitting processing on the distribution map of the phase rotation amount caused by static magnetic field nonuniformity in the standard region; (c) calculating an evaluation value from phase values of a first region that is adjacent to the standard region and phase values of the standard region acquired in the step (b), (d) performing the unwrapping processing on the first region; (e) performing the fitting processing on the standard region and the first region; (f) calculating an evaluation value from the phase values of the standard region and the first region and the phase values of the standard region and the first region acquired in the step (e); (g) giving a desired offset to a nth region; (h) performing the fitting processing on the standard region and the first region; (i) calculating an evaluation value from the phase values of the standard region and the first region and the phase values of the standard region and the first region acquired in the step (h); (j) repeating a desired number of times the steps (g) to (i) while varying the offset value; (k) comparing a plurality of evaluation functions obtained in repetition of the steps (c), (f), (i), and (j); (l) selecting a desired evaluation value on the basis of comparison result in the step (k); and (m) adjusting the phase values of the standard region and the first region on the basis of the evaluation values selected in the step (l).
 9. An image processing method according to claim 8, further comprising a step, in the case where there exist three or more island-like regions, of (n) setting a new standard region that is constituted of the standard region and the first region, the phase values of which have been adjusted in the step (m), and repeating the steps (a) to (n).
 10. A magnetic resonance imaging apparatus according to claim 1, wherein said control means creates an image including at least two separated portions, calculates an approximate function of the phase values of a first region, applies phase values of a second region facing the first region to the approximate function, and thus obtains a water/fat separated image.
 11. A magnetic resonance imaging apparatus according to claim 10, wherein said control means corrects phase rotation caused by the static magnetic field inhomogeneity in the first region, calculates the approximate function, and applies the phase values of the second region to the approximate function.
 12. A magnetic resonance imaging apparatus according to claim 10, wherein the first region is located near the center of the image.
 13. A magnetic resonance imaging apparatus according to claim 10, wherein the first region has the largest area in the image.
 14. A magnetic resonance imaging apparatus according to claim 10, wherein the second region selected from among all regions is the region having a narrow gap with the first region.
 15. A magnetic resonance imaging apparatus according to claim 10, wherein, when phase values of an nth (n is an integer, which is three or more) region are calculated, a function approximating phase values of the first region to the n-1th region are calculated, and phase values of the nth region are applied to the approximate function of the phase values of the first to the n-1th region.
 16. A magnetic resonance imaging apparatus comprising control means for acquiring image data from the object to be examined and obtaining an image by calculation, wherein said control means creates an image including a region that is separated into at least two or more portions, calculates evaluation functions of the respective phase values of a first region and a combined region including a first region and a second region facing the first region, applies the phase values of the second region to the evaluation function from the phase value of the first region, and thus a water/fat separated image of the first region and the second region is obtained.
 17. A magnetic resonance imaging apparatus according to claim 16, wherein said control means corrects phase rotation caused by the static magnetic field inhomogeneity between the first region and the second region, and calculates evaluation functions of the phase values of the respective regions.
 18. A magnetic resonance imaging apparatus according to claim 16, wherein the first region is located near the center of the image.
 19. A magnetic resonance imaging apparatus according to claim 16, wherein the first region has a largest area in the image.
 20. A magnetic resonance imaging apparatus according to claim 16, wherein the second region selected from among all regions is the region having a narrow air gap with the first region.
 21. A magnetic resonance imaging apparatus according to claim 16, wherein, when phase values of an nth (n is an integer,—which is three or more) region are calculated, the phase value of the nth region is applied to the approximating function of the combination of all the regions from the first region to the nth region.
 22. An image processing method for image-processing on a plurality of images so as to obtain an image where water and fat are separated, comprising steps of: (a) performing the unwrapping processing on a first region in an image including a plurality of separated regions for calculating a distribution map of a phase rotation amount caused by the static magnetic field inhomogeneity; (b) performing fitting processing on the distribution map of the phase rotation amount caused by the static magnetic field inhomogeneity in the first region; and (c) calculating an evaluation value from phase values of a second region adjacent to the first region and the phase values of the first region acquired in the step (b).
 23. An image processing method according to claim 22, further comprising steps of: (d) performing the unwrapping processing on the second region; (e) performing the fitting processing on the first region and the second region; and (f) calculating an evaluation value from the phase values of the first region and the second region and phase values of the first region and the second region acquired in the step (e).
 24. An image processing method according to claim 23, further comprising a step of: (g) giving an offset of ±2mπ (m is a natual number) to the phase values of the second region, and thus calculating an evaluation value.
 25. An image processing method according to claim 24, further comprising a step of: (j) comparing evaluation functions of a plurality of evaluation values acquired in the steps (c), (f), and (g), and selecting a desired evaluation value.
 26. An image processing method according to claim 25, further comprising a step of: (h) performing the fitting processing on the phase values of the first region and the second region and phase values of the xth (x is an integer, which is three or more) region. 